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ABSTRACT 

We study the diffe rences and simil arities in the luminosities of bound, infalling and the so- 
called backsplash dGill et al .1120051) galaxies of the Milky Way and M31 using a hydrody- 
namical simulation performed within the Constrained Local UniversE Simulation (CLUES) 
project. The simulation models the formation of the Local Group within a self-consistent cos- 
mological framework. We find that even though backsplash galaxies passed through the virial 
radius of their host halo and hence may have lost a (significant) fraction of their mass, their 
stellar populations are hardly affected. This leaves us with comparable luminosity functions 
for infalling and backsplash galaxies and hence little hope to decipher their past (and differ- 
ent) formation and evolutionary histories by luminosity measurements alone. Nevertheless, 
due to the tidal stripping of dark matter we find that the mass-to-light ratios have changed 
when comparing the various populations against each other: they are highest for the infalling 
galaxies and lowest for the bound satellites with the backsplash galaxies in-between. 
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1 INTRODUCTION 



Ever since iKlypin et all d 19991) and lMoore et al. I GUI) pointed out 
that dark matter simulations of cosmic structure formation lead to 
an excess of subhaloes as compared to the number of observed 
(luminous) satellite galaxies visibly surrounding the Milky Way 
(MW) and M31, the industry for simulating and studying substruc- 
ture in cosmological (dark matter) haloes has boomed. The tension 
has been marginally loosened with the discovery of a substantial 
number of new ultra-faint satellites galaxies in the L ocal Group 
thanks to the SDSS data dAdelman-McCarthvl [2007h : within the 
past couple of years the number of known MW and M31 satel- 
lites has nearly doubled. And taking into account the detection 
limits and the sky coverage of the SDSS survey we will most cer- 
tainly stumble across even more galactic satellites in the near future 
when, for instance, upcoming data from GAIA and panSTARRS 
have been analysed. 

As noted by several groups before dMoore et al.l 120041 : 
iGill et al1l2005l) . there exists a prominent population of galaxies 
that are found outside the virial region of their host at the present 
day, but whose orbits took them inside the virial radius at earlier 
times. While their studies were based upon cosmological simula- 
tions of galaxy clusters, the existence of this "bac ksplash popu- 
lation" has also been fou nd for MW-type objects dWarnick et al .1 
l2008l : lLudlowetal.ll2009h . This raises the question whether or not 



(and how) one can distinguish i nfalling and backsplash galaxies 
from each other. I Gill etaiTd2005h suggested to use the line-of-sight 
velocity distribution: as shown in their Fig.8 the distribution of line- 
of-sight velocities of subhaloes relative to the host (and convolved 
with the 2dF velocity uncertainty of 100 km/sec) is different for 
the infalling and the backsplash population. However, there may 
be a simpler way that does not involve spectroscopy. Since back- 
splash satellites had, at one point in their orbit, a closer approach 
to the central galaxy than infalling satellites, the tidal influence of 
the host must have been stronger for the backsplash population 
then for infalling satellites. Could this difference in tidal forces ef- 
fect the initial stellar population (if existent), and can it be used to 
discrimin a te bet ween the two populations? It has been shown by 
IGill et alj d2005t) that backsplash galaxies loose on average 40% 
of their initial mass when grazing their host. But what about the 
stellar content? As baryons are able to cluster more strongly in the 
centre of the potential well the stars are also more centrally concen- 
trated. Therefore, will the cold baryonic component be safe from 
tidal stripping when the backsplash galaxy (briefly) flies through 
its host? This question is the major motivation for this work. We 
address the issue of separating the three types of galaxies (bound 
satellites, backsplash and infalling) by means of luminosity (and 
possibly mass) measurements only. 
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2 THE SIMULATIONS 

In this Section we describe the simulations used throughout this 
study and the methodology employed to identify host haloes and 
their substructure. 

2.1 Constrained Simulations of the Local Group 

We use the same set of simulations alr eady presented in 
Libcskind ct al. and lKnebe et al. and refer the reader 

to these papers for a more exhaustive discussion and presentation 
of these constrained simulations of the Local Group that form part 
of the CLUES projectQ However, we briefly summarize their main 
properties here for clarity. 

We choose to run our simulations using standar d ACDM ini- 
tial c onditions, that assume a WMAP3 cosmology dSpergel et al .1 
120071) . i.e. a m = 0.24, Cl b = 0.042, fi A = 0.76. We use a nor- 
malization of <r 8 = 0.73 and a n = 0.95 slope of t he power 
spect rum. We used the PMTree-SPH MPI code GADGET2 dSpringell 
l2005h to simulate the evolution of a cosmological box with side 
length of Lbox = 64/r'Mpc. Within this box we identified (in 
a lower-resolution run utilizing 1024 3 particles) the position of a 
model local group tha t closely resembles the real Local Group (cf. 
iLibeskind et al.l20ld) . This Local Group has then been re-sampled 
with 64 times higher mass resolution in a region of 2/r'Mpc about 
its centre giving a nominal resolution equivalent to 4096 3 particles 
giving a mass resolution of otdm = 2.1 x 10 5 /3~'M o for the dark 
matter and m gas = 4.42 x lO 4 /i~ 1 M for the gas partic les. For more 
details we refer to the reader to lGottlober et al.l d2010t) . 

For this particular study we focus on the gas dynamical SPH 
simulation, in which we follow the fe edback and star forma- 
tion rules o f ISpringel & Hernquistl J2003b : the interstellar medium 
(ISM) is modeled as a two phase medium composed of hot am- 
bient gas and cold gas clouds in pressure equilibrium. The ther- 
modynamic properties of the gas are computed in the presence of 
a uniform but evolving ultra-violet cosmic back ground generated 
from QSOs and AGNs and switched on at z = 6 l lHaardt & Madaul 
[19961) . Cooling rates are calculated from a mixture of a primor- 
dial plasma composition. No metal dependent cooling is assumed, 
although the gas is metal enriched due to supernovae explosions. 
Molecular cooling below 10 4 K is also ignored. Cold gas cloud for- 
mation by thermal instability, star formation, the evaporation of gas 
clouds, and the heating of ambient gas by supernova driven winds 
are assumed to all occur simultaneously. 

2.2 The (Sub-)Halo Finding 

In order to identify halos and subhaloes in our simulation we have 
ru n the MPI+OpenMP hybrid halo finder AHlj^] described in detail 
in lKnollmann & K nebe (2009]). AHF is an improvement of the MHF 
halo finder JCrill et al j |2004!) . which locates local over-densities in 
an adaptively smoothed density field as prospective halo centres. 
The local potential minima are computed for each of these density 
peaks and the gravitationally bound particles are determined. Only 
peaks with at least 20 bound particles are considered as haloes and 
retained for further analysis (even though we place a tighter con- 
straint on the number of particles for the present analysis, cf. be- 
low). We like to stress that our halo finding algorithm automatically 

1 http://www.clues-project.org 

2 AMIGA halo finder, to be downloaded freely from 
http: //www. popia. ft.uam.es/AMIGA 



identifies haloes, sub-haloes, sub-subhaloes, etc. For more details 
on the mode of operation and actual functio nality we though refer 
the re ader to the code description paper by iKnollmann & Knebd 
d2009l) . 

Subhaloes are defined as haloes which lie within the virial 
region of a more massive halo, the so-called host halo. We build 
merger trees by cross-correlating haloes in consecutive simulation 
outputs. For this purpose, we use a tool that comes with the AHF 
package called Merger/Tree, that follows each halo (either host or 
subhalo) identified at redshift z = backwards in time. The direct 
progenitor at the previous redshift is the object that shares the most 
particles with the present halo and is closest to it in mass. Again, 
for mo re elaborate details we point to the reader to lLibeskind et al .1 
J2010h . 

2.3 Lighting up Subhaloes 

The stellar populati on synthesis model STARDUST (see 
iDevriendt et al. I 1 19991 , and references therein for a detailed 
description) has been used to derive luminosities from the stars 
formed in our simulation. This model computes the spectral energy 
distribution from the far-UV to the radio, for an instantaneous star- 
burst of a given mass, age and metalicity. The stellar contribution 
to the to tal flux is calcu lated assuming a Kennicutt initial mass 
function dKennicuttll 19981) . 



3 RESULTS 

The prime target of this study is to find possible differences in 
the properties of backsplash, bound and infalling galaxies with re- 
spects to luminosity. We explicitly use the term "galaxies" as we 
focus solely on subhaloes with a luminous stellar component; all 
other objects will be neglected for this particular investigation. We 
further only consider satellites of the (simulated) MW and An- 
dromeda (M31) galaxy; the subhaloes of both these host haloes 
will be stacked in the subsequent plots presented here. In addition 
to the requirement for subhaloes to contain stars we also apply a 
lower mass cut of M > 2 x 10 7 /r'M o roughly corresponding to 
100 particles in total (note that particles have different masses as 
they represent dark matter, gas and stars). 

3.1 The Existence of Backsplash Galaxies 

Before examining the properties of backsplash galaxies we wish 
to first confirm their existence. To this extent we plot in FigQ]the 
closest approach (normalized to the virial radius of the satellite's 
host at the time of minimum distance) vs. its present-day distance 
(normalized to its host's virial radius). The number of objects in 
the respective population are given in the legend. Note that we only 
plot those subhaloes that contain a stellar component. This figure 
contains three distinct parts defining the three different populations. 
First, those subhaloes whose minimum distance equals its present- 
day distance are the infalling population: they are continuously 
falling towards their host. Second, galaxies that entered the virial 
radius of their host and remained inside ever since. Even though 
the host radius is increasing in size since the time a subhalo en- 
tered, we nevertheless find that there are no subhaloes above the 
1:1 line; we therefore conclude that the increase in host radius as 
measured by R'^/R^ is smaller than the ratio D naw /D mm . This 
comes as no surprise as we do not expect satellites to orbit on cir- 
cular orbits (D ]low = D m i„); subhaloes may have (highly) eccentric 
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Figure 2. The orbits of all considered subhaloes. The left panel shows the backsplash galaxies, the middle panel the bound and the right panel the infalling 
population. 
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Figure 1. Minimum distance D m in as a function of present-day distance 
Aiow both normalized to the virial radius of the host at the respective time. 



or bits taking them close to the centre of their host (cf. Figs. 7 and 8 
in iGill et al.l2004h . Third and last, there are galaxies that once were 
inside their host's virial radius but are presently found outside, i.e. 
the backsplash population. 

Fig[TJindicates that we might expect to find of order 40% to be 
backsplash galaxies in the vicinity of the Milky Way and/or M31 
- a percen tage in agreement with previous stud ies of this class of 
objects(cf. IGill et al.l2o"ol ; IWarnick et al.l2008l) . The question now 
is whether or not we will be able to distinguish these populations 
and find the backsplash galaxies, respectively, by quantifying their 
luminosities. 

Before proceeding we would like to add a cautionary remark 
clarifying our terminology: we call subhaloes that are inside their 
host's virial at z = 0, "bound". Those subhaloes that were once 
inside their hosts virial radius but are found at z = outside are 
termed "backsplash". As we can see from Fig[2]this classification 
strongly depends on the redshift used to define the populations. We 
can see that a fair fraction of today's bound population had been 
backsplashed in the past, while probably all of today's backsplash 
galaxies will return and re-enter their hosts at some future time. 
Therefore, the expression "bound" should not be taken literally (in 



terms of energy arguments) but rather as a reference to satellites 
under a prolonged influence of their host while "backsplash" refers 
to satellites under brief influence of their host. Further, please note 
that we require objects to exist both at redshift z = 1.5 and today to 
be part of our sample; there are also subhaloes that were present at 
high redshift but got tidally disrupted and hence did not survive. 

In preparation of the investigation the luminosities in Sec- 
tion !3.2l and baryon content in Section [3~3l we wish to find the time 
where all three populations were still infalling so as to verify the 
correctness of our tracking scheme for subhaloes. To this extent we 
present in Fig[2]the distance from the center of their hosts of all 
backsplash (left panel), bound (middle panel), and infalling (right 
panel) galaxies found and identified at redshift z = 0. The orbits 
have been normalized to the virial radius of the respective host (at 
redshift z) of each galaxy and hence the solid line Am/^hosi = 1 
marks the "entry" (and "exit") point of the satellite into and out of 
the host. While this figure succinctly demonstrate that backsplash 
galaxies clearly entered and exited their host (while infalling galax- 
ies have not yet crossed the virial radius) it also allows us to find 
that point in time in our simulation at which all populations were 
still infalling: this can be seen at a redshift z > 1.5. We will return 
to this redshift later as we expect the properties of galaxies to be 
drawn from the same statistical distribution at that time: no galaxy 
has yet entered the host (or left again) which may (or may not) have 
caused changes in the internal properties and - in particular - the 
luminosities. 



3.2 The Luminosities 

In this section we look at the luminosity of the stellar compo- 
nents of galaxies identified as bound, backsplash and infalling. 
We start by comparing, in Fig[3] the Johnson V-Band luminos- 
ity of the bound satellites to the backsplash and infalling popu- 
lation of galax i es as well as the observational data as taken from 
iKoposov et all d2008h and iMaccio etai] fcOlfjl) . respectively (thin 
solid line, referred to as "Maccio' sample" from now on): these 
data are a combin ation of the volume corrected MW satellite lu- 
minosity func^nJl&no^ovet_al. 2008) augmented with informa- 
tion from lMateol dl998l) and lMaccio et al.l J2010h kindly provided to 
us by Andrea Maccio (personal communication). And even though 
our bound luminosity function agrees with the Maccio data rather 
well, we stress that we included the observational data merely as 
a reference to guide the eye. It is not our prime objective to repro- 
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Figure 3. The luminosity function of subhaloes in the Johnson V-Band. The 
"Maccio" observational data (thin s olid line) is a combinat ion of the volume 

Koposov et all fe008l) augmented with 



corrected MW l uminosity tunc tion IKopos ov et al l 1ZUU8I) augmented with 
information from lMatecl 1 19981) and MacdcTeTalTHoiOl) under the assump- 
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tion of an NFW-like radial distributions of satellites. Note that the compar- 
ison to the observational data is not the prime target of this study and only 
serves as a reference, respectively. 
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Table 1. Kolmogorov-Smirnoff (KS) probabilities p for various compar- 
isons of the luminosity functions presented in Fig[3] 



duce the MW and/or M31 luminosity function of satellite galaxies 
with our simulation. However, a close match (as seen in Fig[3j re- 
assures us that our simulation is not too far fetched and that our 
method for lighting up subhaloes (cf. Section |23t yields credible 
results. The central theme of this paper is the comparison between 
the (numerically obtained) infalling and backsplash population and 
possibilities to decipher them photometrically. 

To better quantify the differences and similarities between 
the respective simulated luminosity functions we applied the 
Kolmogorov-Smirnoff (KS) test that provides us with the signifi- 
cance level p that the null hypothesis that two data sets are drawn 
from the same parent distribution; small values of p 6 [0, 1] show 
that the two cumulative distribution functions (i.e. in our case two 
luminosity functions) are significantly different^] We find that the 
KS probability that the backsplash and infalling distributions have 
been drawn from the same parent function is 67%. The significance 

• We utilized the routine kstwoO as described in IPress et alj il992[) . 



Figure 4. The luminosity function of subhaloes in the Johnson V-Band at 
redshift z = 1.5 (i.e. the redshift at which none of the galaxies has yet 
entered their respective host). 



level is just 11% when comparing the backsplash with the bound 
population and 19% when comparing the infalling with the bound 
satellites. These numbers and probabilities, respectively, have been 
summarized in Table Q] 

In addition to calculating the KS probability p that these distri- 
butions stem from the same parent distribution we also performed 
the experiment of randomly drawing AW-k galaxies from the in- 
falling and bound sample where Mack is the number of backsplash 
galaxies. Comparing the resulting down-sampled luminosity func- 
tions again using a KS test we find a probability p of p « 0.66 
when comparing the backsplash population to the infalling one and 
p « 0,12 when comparing the backsplash to the infalling or bound 
satellites population. All this hints at similarities between back- 
splash and infalling satellites whereas the bound population has 
likely evolved differently. 

Since all bound and backsplash galaxies themselves were, at 
some stage, infalling satellites the differences between the bound 
and backsplash/infalling luminosity function at redshift z = 
should (at least) be lessened when moving to a time where none 
of the objects had entered their host, i.e. redshift z = 1.5 (cf. Fig|2] 
in Section[3j}; the three (cumulative) distributions of luminosities 
at redshift z = 1,5 are presented in Fig|4] Performing the same ex- 
ercises of comparing the various distributions using the KS statistic 
(cf. TableQ]again), we obtain a marginally larger probability for the 
infalling population to agree with the bound and backsplash sub- 
haloes. However, the compatibility between the bound and back- 
splash is actually lowered. In that regards we need to stress that 
the number of satellites - despite combining MW and M31 - is 
not very large (especially not for the backsplash population) and 
hence any (extensive) statistical analysis has to be taken with care. 
Therefore, the probabilities presented here are more indicative of 
possible trends rather than providing hard evidence for similarities 
and/or differences. 

Even though our primary motivation is to find a way to dis- 
tinguish backsplash from infalling satellites that only utilizes pho- 
tometry, we nevertheless present another (observable) correlation: 
the luminosity vs. the velocity dispersion <x„, in Fig[5] As already 
alluded to above when discussing the luminosity function, we also 
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Figure 5. The relation between Johnson V-Band luminosity and subhalo 
velocity dispersion at redshift z = 0. The observational data is taken form 
IWalker et alj <2009l) . 



added observational data (taken from lWalker et alj J20091 their Ta- 
ble 1)) simply to guide the eye. While we also recover the ob- 
served trend in our numerical data, the focus should lie with the 
inf ailing and backsplash galaxies. To better quantify the correla- 
tions between <x„ and My we calculated the Spearman rank coef- 
ficients Q for the observational data it amounts to R s = 0.687 
whereas there appears to be a marginally stronger correlation for 
our bound satellites of R$ = 0.824. However, this "discrepancy" 
is likely due to the different magnitude limits of both the observa- 
tional and numerical data, i.e. the two data sets do only cover the 
same magnitudes in the range M v e [-13, -5]. 

The respective correlation coefficients R s for the backsplash 
and infalling populations are R s = 0.373 and R s = 0.573 respec- 
tively. While there are differences in the strengths of the correla- 
tion we find it difficult to utilize this interdependence to separate 
backsplash from infalling satellites: while the Spearman rank sig- 
nificances Ss are very close to zero for the bound and observational 
data (indicating a reliable determination of the respective R s value) 
they are of order 0.2 for the backsplash and infalling population, 
probably due to the small statistical sample. 

Above, we showed that while all three populations do follow 
the same trend for the My - <x„ relation (coinciding with the trend 
found in observational data), there are nevertheless subtle differ- 
ences in the strength of this correlation, especially for the back- 
splash and infalling population (cf. the different Spearman rank co- 
efficients Rs ). However, the most prominent and well pronounced 
difference can be found when studying the mass-to-light ratios 
M/Ly presented in Fig[6]as a function of V-band magnitude My. 
We stress that the mass M used in this plot is actually the mass 
within the visible radius of the subhalo; we found the distance of the 



4 The Spearman rank coefficient Rs is a non-parametric measure of cor- 
relation: it assesses how well an arbitrary monotonic function describes 
the relationship between two variables, without making any other assump- 
tions about the particular nature of the relationship between the variables 
iKendall & Gibbonslll990l) . Its significance Ss is a value between and 1 
and a small value indicates a significant correlation. We use the IDL routine 
R.C0RRELATEO to calculate both these numbers. 



Figure 6. Mass-to-light ratios (in terms of solar values) as a function of 
V-band luminosity My. The thin solid line represents the observational re- 
lation as found by iMateol j 19981) . The other lines are the best-fit curves 
(with the amplitude as a free parameter) to the bound (dashed), backsplash 
(dotted), and infalling (dot-dashed) population, respectively, with the leg- 
end listing the respective value of the amplitude, too. Note that we used the 
"mass M inside the visible radius" as described in the text for this plot. 



farthe st stellar particle and used the total mass interior to this radius 
as M. IWadepuhl & Springell d2010h already noted that a (substan- 
tial) shift (i.e. A * 5.2) of the observationally determined analytical 
relation 



M/L 



A 2.5 + 



10 7 



(1) 

(M/L) \ L/L ' 

is required jMateo|[l998L A = 1 in there), which is confirmed by 
our data: leaving A as a free parameter and using only the bound, 
backsplash, and infalling satellites we find A = 7.5 ± 0.7 (bound), 
A = 1 1.1 ± 1.2 (backsplash), and A = 19.1+2.4 (infalling) These 
different amplitudes are naturally explained by the differing histo- 
ries and (strengths of) interactions with the host. We will see be- 
low in Section [331 that bound galaxies lost the largest amount of 
their dark matter when compared with the other two populations; 
infalling satellites are in fact still gaining mass through accretion. 
Therefore, taken together with the fact that their luminosities are 
nevertheless still similar, we may infer that the mass-to-light ratios 
should be significantly different. This notion opens up the possi- 
bility to use the relation presented in Fig[6]to separate the three 
populations from each other. In practice this requires not only pho- 
tometric measurements but also a (proxy for the) mass estimation. 

However, using the mass inside the stellar radius also may ex- 
plain the differences found in Fig[6] stars in real satellites may be 
more compact relative to the dark matter than in our simulation, 
and might therefore be less susceptible to tidal stripping (together 
with the dark matter inside the "visible" radius). This would also 
suggest that the differences between our three different populations 
might be smaller if the luminous parts of the satellites were more 
compact. 

5 The fitting to Eq.jT}, i.e. a function M/L(L), has been done using IDL's 
CURVEFIT routine using equal weights for the data points M/L vs. My; the 
reported standard deviations had been returned by CURVEFIT, too. 
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Figure 7. Difference between Johnson V-Band luminosity at present day 
and redshift z = 1.5 as a function of present day halo mass M. 



Figure 8. Ratio of stellar to gas mass at redshift z = 1 .5 as a function of 
present day halo mass M. 



We further like to mention that we not only used the mass in- 
side the stellar radius as a measure for the mass entering the mass- 
to-light ratio. We also applied various other definitions, e.g. the to- 
tal mass inside the virial radius as well as the mass as determined 
from the velocity dispersion under the assumption of virialisation 
and an NFW density profile (both at the virial radius and at 15% 
of the virial radius). While the amplitudes A are certainly differ- 
ent when using different mass estimates, the general trend remains 
unaltered: the M/L ratios for the infalling satellites are shifted up- 
wards with respects to the backsplash population which itself has 
higher ratios than the bound subhaloes. 

However, we also need to bear in mind one of the subtleties of 
halo finding, especially subhalo finding: the definition of mass and 
the edge of a subhalo, respectively. While it is straight forward to 
define an outer edge for an (isolated) field halo (usually defined as 
the radius at which the mean interior density drops below 200 times 
the critical density), the situation is more tricky for subhaloes: they 
have to be truncated at the point where their density profile starts 
to rise again due to the host's background density. Therefore, the 
same subhalo placed inside and outside of a host will have differ- 
ent masses due to the nature of (our) halo finding technique. This 
explains at least in part the offset in the mass-to-light ratios for 
bound/backsplash and infalling galaxies: the infalling ones have 
in general higher masses. And part of the gap between bound and 
backsplash may be explained by the same phenomenon, though not 
all of it; there certainly is no uncertainty that bound galaxies have 
lost more mass than backsplash subhaloes. 

The differences between the luminosities of the populations 
at redshift z = and the (marginally) more pronounced similari- 
ties at a time where all populations were still infalling (i.e. redshift 
z = 1.5) calls for a closer look at the evolution of satellite galaxy 
luminosity. To this extent we plot in Fig|7]the difference between 
the Johnson V-Band luminosity My at redshift z = and at red- 
shift z = 1.5 as a function of the total bound halo mass M for each 
galaxy considered in this paper, using different symbols for the dif- 
ferent populations (stars for backsplash, plus-signs for bound, and 
diamonds for infalling galaxies). 

Fig|7]now shows several things. For a substantial number of 
satellites (especially the backsplash and infalling population) we 



only observe a "constant" decrease in luminosity of approx. 1.5 
magnitudes. However, the luminosity of the bound galaxies drops 
significantly - especially on the low-mass end - while some of the 
higher mass ones gain luminosity. As luminosity is directly linked 
to stellar content we are left with the question of how these differ- 
ences relate to changes in the stellar population and/or removal (or 
gain) of star particles from a subhalo. We study these issues in the 
following subsection. 

Studying luminosities is also closely related to colours, i.e. 
ratios of luminosities in different wave-bands. It therefore appears 
natural to ask the question - and use the data available to us - to 
have a closer look at differences in colours for our three popula- 
tions. When plotting the B-V colour as a function of halo mass M 
(not presented here) we observe that there are practically no differ- 
ences at all amongst the various subhaloes and populations, respec- 
tively. Neither is there are correlation with mass. Colour appears 
unaffected when categorizing galaxies as bound, infalling or back- 
splash. 

3.3 The Baryon Content 

Before investigating the stellar component directly we would like 
to start with a few words on the subhalo gas content. We find that 
hardly any of the subhaloes under consideration contain a signif- 
icant gas content at redshift z = 0. When expressed in terms of 
the cosmic baryon fraction the amount of mass in gas is of order 
< 10~ 4 for more than 90% of the subhaloes. However, their stel- 
lar mass fractions (again in terms of the cosmic baryon fraction) is 
> 10~ 4 for all of them (which is a direct outcome of restricting our- 
selves to a sample of subhaloes that contain a stellar component). 
The situation though is rather different at redshift z = 1.5 where we 
find that all of the progenitors of the subhaloes not only contained 
gas but the fraction of mass in gas is on average a factor two higher 
than that in stars which can be verified in Fig[8] Note that in this 
figure we plot the present day mass on the x-axis and the ratio of 
stellar-to-gas mass at redshift z = 1.5 on the y-axis. The fact that 
none of the gas is left at redshift z = indicates that either the gas 
has been converted into stars through the process of star formation, 
or the gas has been stripped/removed through interactions with the 
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Figure 9. Ratio of stellar mass at present day and redshift z = 1.5 as a 
function of present day halo mass M. 



Figure 10. Ratio of total bound mass at present day and at redshift z = 1 .5 
as a function of present day halo mass M. 



host halo (e.g. ram pressure stripping) or other influences prior to 
infall. If the former is true we would expect to observe an increase 
in stellar mass since redshift z = 1.5, unless there is a conspiracy 
at work: the existing stars may be stripped at the same rate as star 
formation may convert gas into new stars, leaving the number of 
stars unchanged. However, this scenario is rather unlikely. We have 
also checked for the influence of the cosmological UV background: 
recall that in our simulations the thermodynamic properties of the 
gas are computed in the presence of a uniform but evolving UV 
cosmic background generated from quasi- stellar objects and activ e 
galactic nuclei and switched on at z = 6 dHaardt & Madau 0*9961) . 
This prescription leads to an evaporation of gas in o bjects below 
a certa in mass threshold M c (z) as given by Eq. (6) in lHoeft et all 
( 120061) . When plotting the mass accretion histories of all our sub- 
haloes under inv estigation here and comparing it to the aforemen- 
tioned formula in lHoeft et all j2006l) . we find that all the backsplash 
and infalling galaxies are in fact below the evaporation limit. For 
the bound objects we find that 1/3 are, at redshift z = 1.5, above 
that mass limit. However, they also drop below it by z = (with the 
odd one remaining above). We therefore conclude that we should 
not be surprised to be left with subhaloes that contain hardly any 
gas at z = 0, due to photo-evaporation by the UV background. 

The question now is, whether or not we find an evolution of 
the stellar component between redshifts z = 1.5 and z = 0. We 
therefore plot the ratio of the stellar content at these redshifts as a 
function of (present-day) subhalo mass in Fig|9] We observe that 
the backsplash (as well as the infalling) subhaloes hardly lost any 
stars since z = 1.5. Note that our simulations do not model stel- 
lar mass loss and hence the stellar mass remains constant when no 
star particles are stripped or newly created. However, this is still in 
agreement with the evolution of the luminosity as found in Fig(7] 
subhaloes with a constant number of stellar particles merely evolve 
passively from z = 1.5 to z = due to stellar ageing. Nevertheless, 
the lower-mass subhaloes of the bound population did loose a sub- 
stantial amount of stars while some of the higher-mass ones gained 
(or formed) stars. Therefore, a picture is now emerging, that while 
gas has been efficiently stripped, the stellar component remained 
more or less unaffected - at least for the infalling and backsplash 
population which are of prime interest in the present study. 



As we expect the stellar component to be concentrated at the 
centre of a subhalo, the previous finding on stellar mass loss for 
bound subhaloes immediately leads to questions regarding the na- 
ture of mass loss in general. To this extent we show in FigllOl 
the ratio of total bound masses M (again as a function of today's 
mass) at redshifts z = and z = 1.5. We note that outside the 
influence of a host halo, subhaloes behave like field haloes and 
grow in mass through accretion processes; this is clearly confirmed 
for the infalling population (albeit with the exception of two ob- 
jects). We also observe mass loss via tidal stripping, especially for 
the bound subhaloes. And while backsplash g alaxies may at ti mes 
loose as much as 40% of their original mass dGill et alj|2005r) 
also find the odd backsplash galaxy in our particular sample that 
gained mass. Nevertheless, the picture is more or less clear: even 
though backsplash galaxies loose mass, their stellar component re- 
mains unaffected. This is not the situation for bound galaxies that 
loose both dark matter and stars due to the tidal interactions with 
the host, as expected. The picture drawn here therefore naturally 
explains the differences in the (amplitude of the) M/L v ratios even 
though for that particular study only the "mass inside the visible 
radius" has been considered: when using the total bound mass (not 
presented here) we recover the same relations amongst the different 
subhalo populations with the ratios in amplitude unchanged; how- 
ever, the absolute value of the amplitudes is more than a factor two 
higher. 



4 DISCUSSION AND CONCLUSIONS 

In this study we set out to examine the differences of the lumi- 
nosities of backsplash, bound, and infalling satellite galaxies in a 
constrained cosmological hydrodynamical simulation of the Local 
Group. Our prime question is: Is it possible to distinguish these 
different population by mere photometry? While we find marginal 
differences in the bound vs. backsplash/infalling galaxies, the two 
populations residing in the outskirts of the host halo appear to have 
strikingly similar properties in terms of luminosity. The time back- 
splash subhaloes spent under the influence of the host is therefore 
not long enough to affect the stellar component: they loose mass, 
but primarily dark matter and/or gas particles are stripped - the star 
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particles remain more or less unaffected by the host's tidal field. 
Therefore, their luminosity function and luminosities in general re- 
main akin to the infalling population. 

Nevertheless, when allowing for not only photometric infor- 
mation but also adding "mass" to our analysis, we found that the 
mass-to-light ratios (as a function of magnitude) are significantly 
higher in infalling than in backsplash galaxies, which are in turn 
both higher than for bound sat ellites. Fitting the observationally 
determined relation presented in lMateol d 19981) for M/L v vs. M v by 
leaving the amplitude as a free parameter we find differences in the 
amplitude of a factor of 1.5 and 2.5 between backsplash and bound 
and infalling and backsplash galaxies, respectively. We note, how- 
ever, that part of this shift can be explained by (our method of) halo 
finding and certain endemic limitations when comparing field and 
subhaloes. The radial extent of a subhalo has to be truncated due to 
the embedding within the host's background and hence has a lower 
mass than in the case when the same subhalo is found in isolation 
(i.e. exterior to a host halo) even though we explicitly only con- 
sidered the "mass inside the visible radius" for this particular part 
of the investigation. We also need to acknowledge that the original 
relation had to be shifted by a factor of 7.2 to bring it into agree- 
ment with our numerical dat^fl which nevertheless is not the prime 
target of the present paper and its explanation left to a future study, 
respectively. 

Even though there still remains a lot to be quantified, we be- 
lieve that this difference may provide a new window on distinguish- 
ing between infalling and backsplash galaxies that could be applied 
to observational data. Its origin is readily explained by the fact that 
while backsplash and bound galaxies both lose mass the mass loss 
is greater for bound than for backsplash galaxies; therefore, if the 
stellar population is unaffected (as found in our simulations) we 
observe an enhanced decrease in the mass-to-light ratios for bound 
galaxies and - more importantly for our purposes - a decrease when 
comparing infalling against backsplash. 

However, the apparent discrepancy between the simulation 
presented and used here and the observational data yet remains 
unexplained. It could be possible that stars in real satellites are 
more compact relative to the dark matter than in our simulation, 
and might therefore be less susceptible to tidal stripping (together 
with the dark matter inside the "visible" radius). But this would 
also suggest that the differences between our three different popu- 
lations might be smaller if the luminous parts of the satellites were 
more compact closing the aforementioned "window" again. In that 
regards, we remind the reader that we not only used the mass in- 
side the stellar radius as a measure for the mass entering the mass- 
to-light ratio. We also applied various other possibilities, e.g. the 
total mass inside the virial radius as well as the mass as determined 
from the velocity dispersion under the assumption of virialisation 
and an NFW density profile (both at the virial radius and at 15% of 
the virial radius). While the ratios of the M/L curves are certainly 
different when using different mass estimates, the forcited trend re- 
mained unaltered. 

A closer inspection of Fig[2] reveals that most of the back- 
splash galaxies fell into their host at approximately the same time. 
When studying the distribution of infall times (not shown here 
though) there appears to be a continuous infall of bound galax- 
ies whereas the backsplash objects all cluster at about redshift z ~ 
0.55. As pointed out by several other authors recently, subhaloes 



6 Note that jWadepuhl & Springel 1 201(f) also required a shift by a factor of 
5.2 for their simulation data. 



ma y have the tendency to fall into (Milky Way like) hosts in groups 



(cf. Klimentowski et al. 


20ld:lLi & Helmil2009l;lD'Onehia & Lake! 


l2008l:lLi&Helmill2008 


). Hence could it be that all our backsplash 



galaxies are part of a larger group? We explicitly checked for this 
conjecture by studying their 3D orbits and cannot confirm it: our 
backsplash galaxies come from various directions yet fall in at a 
similar redshift. However, we also need to acknowledge that these 
directions are not random but rather correl ated - however, t his ha s 
been studied in detail in a companion paper Libes kind et al. I J2O10h . 
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